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ABSTRACT 

o 

04 ' Context. Tidal dissipation in late-type stars is presently poorly understood and the study of planetary systems hosting hot Jupiters can 

^ , provide new observational constraints to test proposed theories. 

. Aims. We focus on systems with F-type main-sequence stars and find that the recently discovered system CoRoT-1 1 is presently the 



best suited for such a kind of investigation. 

Methods. A classic constant tidal lag model is applied to reproduce the evolution of the system from a plausible nearly synchronous 
00 , state on the ZAMS to the present state, thus putting constraints on the average modified tidal quality factor {Q'^ of its F6V star. Initial 

conditions with the stellar rotation period longer than the orbital period of the planet can be excluded on the basis of the presently 
observed state in which the star spins faster than the planet orbit. 

Results. It is found that 4 x 10^ ^ (Q'J i2x 10^, if the system started its evolution on the ZAMS close to synchronization, with an 
uncertainty related to the constant tidal lag hypothesis and the estimated stellar magnetic braking within a factor of » 5-6. For a 
non-synchronous initial state of the system, (QD ^ 4x 10* implies an age younger than ~ 1 Gyr, while (Q'J ^ 2x 10^ may be tested by 
(— I ' comparing the theoretically derived initial orbital and stellar rotation periods with those of a sample of observed systems. Moreover, 

^ ^, we discuss how the present value of Q'^ can be measured by a timing of the mid-epoch and duration of the transits as well as of the 

I . planetary eclipses to be observed in the infrared with an accuracy of ~ 0.5 -Is over a time baseline of ~ 25 yr. 

O ' Conclusions. CoRoT-1 1 is an highly interesting system potentially allowing us a direct measure of the tidal dissipation in an F-type 

star as well as the detection of the precession of the orbital plane of the planet that provides us with an accurate upper limit for the 
C/5 , obliquity of the stellar equator. If the planetary orbit has a significant eccentricity (e ^ 0.05), it will be possible to detect also the 

precession of the line of the apsides and derive information on the Love number of the planet and its tidal quality factor. 

Key words, planetary systems - planet-star interactions - binaries: close - stars: late-type - stars: rotation 

> 

1 . Introduction lar coordinates with the origin in the centre of the star, w/,,, is the 

t — ..J... , ■ . , . tidal frequency in that frame, *P;„, the amplitude of the compo- 

1.1. Tidal dissipation theories ^^^^ ^^^^.^^ ^ azimuthal order m, and f the time. The tidal 

(Sj Tidal dissipation in close binaiy systems with late-type compo- frequency is given by ^ In - mQ, where n = 27r/Porb is the 

^ nents is generally constrained by the ranges of orbital periods mean motion of the binary, Po,-b being its orbital period. Waves 

O ! correspo nding to circular orbits as observed in clusters of differ- are expected to be excited with the different frequencies w,,„ cor- 

^ . ent ages. iQgilvie & Lini d^O?) review recent observations and responding to the various components of the tidal potential and 

■ conclude that the equilibrium tide theory is insufficient to ex- their amplitudes will depend on that of the exciting component 

. !-( plain binary circularization by at least two orders of magnitude, and on the response of the stellar- interior For a nearly circular 

X ' Therefore, in addition to the dissipation of the kinetic energy of orbit, the / = m = 2 component is the dominant one and is re- 

H ■ the flow associated with the tidal bulge, which is considered in sponsible for the synchronization ophe stellar rotation with the 

. ."^ ! ttie equilibrium tide theory (e.g., IZafflfTOTTl [1989), other ef- orbital motion (e.g., |Ogilvie & Lin||2004,). 

fects must be included. The dynamical tide theory treats the dis- The efficiency of tidal dissipation is usually parameterized 

sipation of waves excited by the oscillating tidal potential in the by a dimensionless quality factor Q proportional to the ratio of 

stellar- interior whose kinetic energy is ultimately extracted from the total kinetic energy of the tidal distortion to the energy dissi- 

the orbital motion. For simplicity, we shall assume fliat stars are pated in one tidal period 2n/ib (e.g., Zahn 2008). In the theory, 

rotating rigidly. Q always appears in the combination Q' = (3/2)(2/^2), where 

We consider a reference frame rotating with the stellar an- k is the Love number of the star that measures its density strat- 

gular velocity Q = 27r/P„t, where P„t is the stellar rota- ificatioifl Therefore, the smaller the value of Q', the stronger 

tion period. In that frame, the tidal potential experienced by the tidal dissipation. In general, Q' depends on I, m, and the 

flie star can be written as a sum of rigidly rotating compo- tidal frequency cb, thus a rigorous treatment of the tidal dissi- 

nents proportional to the spherical hai-monics Yi„,{9,4>), viz., pation should consider the sum of the effects associated with the 

Re [*l'/„,r'}','„(0, 0) exp(-/<y/„,f)], where (r,0,(p) are spherical po- 

Note that k2 is twice the apsidal motio n constant of the star, often 
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different tidal components each having its specific Q . In prac- 
tise, we adopt a single value of Q which represents an aver- 
age of the contributions of the different components. Moreover, 
we also average on the tidal frequency, which means averaging 
along the evolution of a given system because the tidal frequency 
decreases with time and goes to zero when tidal dissipation has 

circularized and synchronized the binary. 

The observations reviewed by lOgilvie & Lin I ( l2007h indi- 
cate that an average Q ranging between 5 x 10^ and 2 x 10^ 
is adequate to account for the circularization of late-type main- 
sequence binaries. Such low values require an efficient tidal dis- 
sipation mechanism that Qgilvi e & Lin I ^2007), moving along 
the lines of previous work, propose to be the damping of inertial 
waves in the stellar interior These waves have the Coriolis force 
as their restoring force and are excited provided that the tidal 
frequency w satisfies the relationship; 

|a>| < 2Q. (1) 

The corresponding Q has a remarkable dependence on the tidal 
and stellar rotation frequencies owing to the complex details of 
wave excitation and dissipation that are still poorly understood 
(cf. also lGoodman & Lackner 120091) . The main point, which can 
be regarded as firmly established, is that Q decreases by 2 - 4 
orders of magnitude when |w|/Q < 2 with respect to the case 
when |fa;|/Q > 2, because the excitation of inertial waves is for- 
bidden in the latter case and only the damping of the equilibrium 
tide contributes to the dissipation. 

In view of the present uncertainties in the dynamical tide the- 
ory, a precise determination of the tidal dissipation in binary sys- 
tems with well known parameters is highly desirable. The case 
of F-type main-sequence stars is particularly challenging from 
a theoretical viewpoint because their internal structure consists 
of a thin outer convective zone and a radiative interior hosting 
a small convective core at the centre of the star. Since the prop- 
agation and dissipation of inertial waves are remarkably differ- 
ent in the convective and radiative zones, the study of F-type 
stars provides a critical test for the theory. As the mass of the 
outer convection zone decreases rapidly with increasing stellar 
mass between 1.2 and 1.5 M©, the value of Q is expected to 
increase by 3 - 4 orders of magnitude within this mass range 
(lBarker&OpiIviell2009h . 



7.2. Testing tidal tiieory witli planetary systems 

A new opportunity to test the tidal theory comes from the 
star-planet systems, in particular those containing hot Jupiters. 
Systems with transiting planets have the best determined stel- 
lar and planetary parameters and are p articularly suited to study 
tidal dissipation (see, e.g., ICarone & Patzold1l2007h . F-type 
host stars having a mass M > 1.2-1.5 M© evolve quite 
rapidly during their main-sequence lifetime, thus improving sig- 
nificantly their age estimate from model isochrone fitting in 
comparison with lower mass stars. A good age estimate is im- 
portant to constrain the average value of Q' by modelling the 
tidal evolution of a particular system (cf. Sect. 13.31 1. In Table [1] 
we list the presently known transiting systems with a star having 
an effective temperature Tetf > 6250 K which corresponds to a 
spectral type earlier or equal to F8 V. 

The columns from left to right list the name of the system, 
the effective temperature Teff, the mass M and the radius R of 
the star, its rotation period P^ot, as derived from the observed 
spectroscopic rotation broadening v sin / and the estimated stel- 
lar radius assuming a equator-on view of the star, the orbital pe- 
riod Porb, the mass Mp and radius Rp of the planet, the timescale 



for the synchronization of the stellar rotation T^yn, and the refer- 
ences. To compute the synchronization time, we assume that the 
entire star is synchronized as is customary in tidal theory and 
is suggested by the tidal evolution of close binaries observed in 
stellar clusters of different ages. The synchronization timescale 
is then a measure of the strength of tidal dissipation in the star 
and is computed according to the formula: 



dD. 
dt 



ly^Ql \M j \a) \n) 



GM 



(2) 



where yR ^ 0.22 R is the gyration radius of the star (Siess et aTl 
2000), 2s its modified tidal quality factor, here assumed to be 
Q's - 10^, a th e semimajor axis of the orbit, and G the gravitation 
constant (see iMardUng &Linll2002h . Equation (H) is valid for 
circular orbits and when the spin axis is aligned with the orbital 
angular momentum. In this regard, the values given here must 
be considered as estimations and for illustration purpose only, as 
high eccentricity and/or obliquity have been measured for some 
of these systems. 

According to the dynamical tide theory by lOgilvie & Lin I 
(I2OO7I) . the stars experiencing the strongest tidal interaction are 
those with n/Q = PmtlPorh 5 2 because they have < 2 

when th e I — m - 2 cornponen t of the tidal potential is con- 
sidered ("Barker & Qgilvie II2009I) . For such stars, the orbital de- 
cay (or expansion) due to the tidal interaction can in princi- 
ple be obs erved. Among those systems, the best candidate is 
CoRoT-11 (" Gandolfi etaLll2010l) because it has a tidal synchro- 
nization timescale in between those derived for systems like 
HD 197286/WASP-7 or HAT-P-09, i.e., much longer than the 
main-sequence lifetime of the system, and those of, e.g., OGLE- 
TR-9 or WASP-3 that are shorter than the expected ages of the 
systems, implying that the synchronous final state has possibly 
already been reached. Other systems, e.g., HD15082/WASP-33, 
have a star so massive that a Q'^ too large to be measurable is ex- 
pected, or show a remarkable misalignment between the stellar 
spin and the orbital angular momentum that makes the deriva- 
tion of the rotation period from the v sin / quite uncertain, as in 
the case of XO-4. 

In view of the peculiar characteristics of the CoRoT-1 1 sys- 
tem, we shall consider it for a detailed study of the tidal evolu- 
tion. We shall derive constraints on the average Q'^ value of its 
F6V star by considering a possible initial state for the system 
when its star settled on the zero-age main-sequence (hereafter 
ZAMS; see Sect. 12.4b . Moreover, we shall demonstrate how the 
present value of Q'^ can be directly measured with suitable transit 
observations extended on a time interval of a few decades. 



2. Method: tidal evolution model and initial 
conditions 

2.1. The constant quality factor approximation 

We adopt the classic equiUb rium tidal model after iHut I (Il981h . 
in the formulation given by iLeconte et al. I (l2010l) . It assumes 
that the energy of the tidal perturbation raised by the tidal po- 
tential is dissipated by viscous effects producing a constant time 
lag At between the maximum of the tidal potential and the tidal 
bulge. When Q» I, 2"' - ^^^t ~ 26(0)), where 5(w) is the lag 
angle between the maximum of the deforming potential and the 
tidal bulge. Since At is assumed to be constant during the tidal 
evolution, Q and Q' vary as w"' . To model the evolution of a sys- 
tem by assuming a constant (or, bette r, an average) Q', we need 
to make an approximation. Following i Leconte et al. (1201 Ol) . we 
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Name 


eff 




M 




K 




^rot 






(K) 




(Mo) 




/D \ 




(days) 




CoKol-l 1 


ii/i /in 1 
044U ± 


1 on 
izU 


i.Z/ ± 


n nc 
0.05 


1.3/ ± 


n nQ 
O.Oj 


1 1 
I. /J ± 


n OA 
0.z6 


I I *rp Yf p./: 

HAl-r-Uo 


03 /U ± 


on 
80 


1 on 1 
i.Z9 ± 


n nA 
O.OD 


1.4d ± 


n nA 
O.Oo 


o /in 1 
8.49 ± 


1 lA 

1.34 


II AT D r\n 
HAl-r-U/ 


oJSU ± 


on 
80 


^ An 1 
i.4/ ± 


n m 
0.0/ 


1 Oyl 1 

l.o4 ± 


n 1 T 
U. 1 / 


0/1 CI 

z4. J 1 


1 c cn 
t 5.59 


rlAl-r-UV 


ojjU ± 


1 


1.Z8 ± 


nil 


1 1^ _i_ 

I.JZ ± 


n m 
U.U/ 


C A 1 _i_ 

J. 61 ± 


n TO 
U. /6 


IT AT D 1 /I 

rlAl-r-14 


OOUU ± 


nn 


1 "in _i_ 


n f\i 


1 AH _i_ 

1.4/ ± 


n n^ 
U.U3 


O OC _i_ 
6.6J ± 


n oc 


U AT D 1/1 

rlAl-r-Z4 


bi 1 i ± 


on 


1 1 n _i_ 


n n/1 
U.U4 


1 _i_ 

l.iz ± 


n m 
U.U/ 


A AT _i_ 

6.6/ ± 


n AO 
U.D6 




oz9U ± 


i iU 


i.iz ± 


n no 


1 /IT 1 

1.4z ± 


n nA 
U.UO 


Q 1 /I 1 
3. 14 ± 


n Qn 


1 1 r> 1 cnoo 
rlJJiDUoZ 




1 nn 


1 ^n _i_ 


n f\i 


] A A ^ 

1.44 ± 


n {\i 
U.UJ 


n o 1 _i_ 
U.61 ± 


n 1 1 
U. 1 1 


riuiy 1 loo 


o4UU ± 


1 nn 
100 


loo 1 


n 1 A 
O.lo 


l.z4 ± 


n nc 
U.Uj 


"3 AO 1 
J. 66 ± 


n cn 
U.D9 


Kepler-IJj 


bly 1 ± 


An 
60 


V.i 1 ± 


n nA 
O.Oo 


1. /y ± 


n nA 
U.Uo 


1 o n 1 
16.91 


1 /ion 
± 4.6U 


Kepler- Uo 


OZlj ± 


1 cn 
150 


1 1 1 


n m 
0.0/ 


1 /in 1 

i.4y ± 


n nA 
U.Uo 


T 1 A 1 

/.lo ± 


n TO 
U. /6 


r^/^T c TO T n 


Arm 1 


An 
60 


i.JZ ± 


n no 
O.Oo 


1 1 

l.jj ± 


n n/1 
U.U4 


1 m 1 
1.9/ ± 


n m 
U.U/ 


\17A CD m 

WAiSr-Uj 


o4l)U ± 


1 nn 
100 


l.z4 ± 


n nn 
0.09 


1.31 ± 


n nn 
U.U9 


/I nc 1 
4.95 ± 


n n 1 
U.91 


\17^A CD IT 

WAc»r-lz 


oZjU ± 


1 cn 
150 


I.j5 ± 


n 1 /I 


1 C7 1 
[.J / ± 


n m 
U.U/ 


TA 1 O 
JO. 1/ 


1 A n nQ 
t 49. Uj 


IT 7 A CD 1/1 


6475 ± 


100 


1.21 ± 


0.12 


1.31 ± 


0.07 


23.68 


t 6.35 


WASP- 15 


6300 ± 


100 


1.18 ± 


0.12 


1.48 + 


0.07 


18.69 


t 13.64 


WASP-17 


6550 ± 


100 


1.20 ± 


0.12 


1.38 + 


0.20 


7.76 ± 


2.49 


WASP- 18 


6400 ± 


100 


1.25 ± 


0.13 


1.22 + 


0.07 


5.60 ± 


1.11 


XO-3 


6429 ± 


100 


1.21 ± 


0.07 


1.38 + 


0.08 


3.81 ± 


0.23 


XO-4 


6397 ± 


70 


1.32 ± 


0.02 


1.56 + 


0.05 


8.97 ± 


0.80 



^orb 

(days) 



(Mj) 



(Rj) 



(Gyr) 



2.994 2.33 ± 0.34 1.43 ± 0.03 8.673 

3.853 1.06 ±0.12 1.33 ±0.06 28.38 

2.205 1.78 ±0.01 1.36 ±0.14 0.166 

3.923 0.78 ±0.09 1.40 ± 0.06 205.8 

4.628 2.20 ±0.04 1.20 ±0.58 15.21 

3.355 0.69 ± 0.03 1.24 ± 0.07 58.49 

5.633 8.62 ± 0.55 0.98 ± 0.04 4.379 

1.220 4.10 ±4.00 1.50 ± 0.05 0.347 

4.955 0.96 ± 0.20 0.93 ± 0.04 553.7 

3.548 2.11 ± 0.06 1.43 ± 0.05 1.009 

3.523 0.60 ±0.19 1.42 ±0.06 61.17 

2.486 4.50 ±1.50 1.61 ±0.04 0.927 

1.847 1.76 ± 0.09 1.29 ± 0.09 1.002 

1.091 1.41 ± 0.09 1.79 ± 0.09 0.012 

2.244 7.34 ± 0.50 1.28 ± 0.08 0.017 

3.752 0.54 ± 0.05 1.43 ± 0.08 22.41 

3.735 0.49 ± 0.06 1.74 ± 0.24 128.7 

0.941 10.30 ±0.69 1.11 ±0.06 0.002 

3.192 11.79 ± 0.59 1.22 ± 0.07 0.814 

4.125 1.78 ± 0.08 1.34 ±0.05 11.32 

Table 1. Parameters of the transiting planetary systems havin g stars with Te.n > 6250 K; Mj - 1 .90 x 10^^ kg and Rj = 
indicate the mass and the radius of Ju piter. Referen ces: A m09 : | Am mler- von Eiff etal.l(l2009l): An09 : [Anderson et al. I 



(2009a); He09b: Hellier at al. (2009b); He09: Hebb et al 



[Bakos et al, (2007 J); CalO:[C ameron et al. I (l2010h; Ga08:lGan dolfi et al. I (|2010l); Gi08:lGibson et al. I(l2008h; He09a 

/onn«_.. ,,.n«,.. ,T_„:.„ /onn^i... T T.n« . 't T ./onn^.. . _ 1 n . U.^..: . _ J, '^2010'); - ' 



(2009); JelO 

[Johns-Krull et al. (2008); Jo09: Johnson et al. (2009); KilO: Kipping et al 
(2008); Mc08 
IPoUacco et al 



Jenkins et al 
(2010); KolO 



Refs. 



GalO 

To08, No08 
Na09 

Am09, Sh09 
TolO, SilO 
KilO 

Ba07, Lo08 

CalO 

He09a 

KolO 

JelO 

Sn09 

Gi08, Po08 
He09 

Jos08, Jo09 
We09 
An09 
He09b 
Joh08, Wi09 
Mc08, NalO 

7.15 X 10^ m 
20ia); Ba07: 



Jos08: Joshi et al 



.Hellier et all 
2008"); Joh08: 
Koch et al. (2010); Lo08: Loeillet et all 
McCullough et al. (2008); Na09: Narita et al. (20091); Nal OilNarita et al. 1(12 010); No08: Noves et al. (2008); Po08: 
(■2008>); Sh09:,Shporer et a l. (2009); SilO: Simpso n et aO (l2010l) ; Sn09nSnellen et al. , (.2009.) ; To08 : .Torres et all 



(I2008h : TolO: [Torres et al. I (120101) ; We09: West et al. (2009); Wi09: IWinn etaL 



(l2009t) 



assume for the star and the planet A;2sAfs = 3/(2«obs2s) 
A:2pAfp = 3/(2nobs2p), where the indexes s and p refer to the 
star and the planet, respectively, and Wobs is the present mean or- 
bital motion of the system corresponding to the orbital period 
f orb = 2.99433 days. As we shall see (cf. Sect. O, the mean 
motion n varies by a factor of s; 4 - 5 during the evolution of 
the system, but we shall neglect such a change because we are 
interested in deriving the order of magnitude of the average 
along the evolution. 

The theory of dynamical tides predicts that the dissipation is 
a complex function of the tidal frequency, the properties of the 
interior of the body, and, possibly, the strength of the tidal po- 
tential because of non-linear effects. All such dependencies are 
lumped together into the coefficients Q'^ and Q'^ and in princi- 
ple could be included in an equilibrium model provided that the 
relationship between their variations along the evolution of the 
system and the time lag value were known. In view of our igno- 
rance of the processes contributing to the tidal dissipation in stars 
and p lanets (cf., e.g., Ogilvie & Lin 2004; Goodman & Lackner 
l2009l) and the presently limited observational constraints, we 
shall consider constant, i.e. average, values of Q'^ and Q'^ all 
along the evolution a nd integrate the equations in Sect. 2.2 of 
iLeconte et al. I (l2010h accordingly. 

2.2. Range of considered Q[ 

The dynamical tide theory of lBarker & Qgilviel (l2009l) gives the 
dependence of Q'^ on the tidal and the rotation frequencies of the 
star when the dissipation of inertial waves is regarded as the main 
source of energy dissipation - see their Fig. 8, a nd consider that 
Q[ PC QT'^ for a fixed ratio {hjO. (cf. Sect. 3.6 of lOgilvie & Lin I 
l2007h . In view of the simplified treatment of the interaction of 



such waves with turbulent convection and radiative-convective 
boundaries, their results can be regarded only as an approxima- 
tion of the complex dependence of Q'^ on the relevant parame- 
ters. For a star of 1.2 and Piot ~ 1.5 days, the average Q'^ for 
|<£i|/Q < 2 can be roughly estimated to be 5 x 10^ with an uncer- 
tainty of at least one order of magnitude and some preference 
for the lower bound because of the approximate treatment of 
the above mentioned processes. Outside such a frequency range, 
2s ~ 10"^, i.e., the tidal dissipation is reduced by at least three 
orders of magnitudes. 



2.3. Range of considered Q'p 

An estimate of the tidal quality factor of the hot Jupiters may be 
based on astrometric observations of the Jupiter-I o system that 
give gp ^ (1.36 + 0.21) X 10^ for Jupiter (see ' Lainev et all 
2009). Nevertheless, the quality factor of hot Jupiters can be dif- 
ferent from that of Jupiter owing to their different internal struc- 
ture. Assuming the observed eccentrit y of the orbits of s evera l 
hot Jupiters to be of primordial origin, lMatsumura et al. I (l2008h 
infer 10^ < gp < 10^, with most of the planets having Q'^ < 10^. 
Therefore, we shall consider 10'' < gp < 10'^ when integrating 
the equation of tidal evolution of our system. The timescale to 
achieve the pseudosynchronization of the planetary rotation in 
the case of CoRoT-11 and of the other systems containing hot 
Jupiters is of the order of 10^ yr when gp ~ 10^ (cf. Sect. 2.2 
of iLeconte et al. 1 l20Toh . Even if Q'^ ~ 10^ such a timescale 

is shorter than 3 x 10^ yr, i.e., much shorter than the typical 
timescales for the tidal evolution of the orbital parameters and 
the stellar spin. Therefore, we shall consider the planet to be al- 
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6000 6200 6400 6600 6800 7000 7200 7400 
Te„ (K) 

Fig. 1. Upper panel: The synchronization parameter n/Q vs. the 
effective temperature of the star Teff in transiting planetary sys- 
tems with Teff > 6000 K. WASP-12 with «/Q ^ 33 has been 
omitted for the sake of clarity; lower panel: an enlargement of 
the lower portion of the upper panel, to better show the domains 
close to n/fi = 1 and n/Q = 2, which are marked by horizontal 
dashed lines. The names of the systems are reported in both pan- 
els, although they are omitted for n/Q, < 4.0 in the upper panel 
to avoid confusion. 



ways in a state of pseudosynchronization which simplifies the 
integration of the tidal equations. 

2.4. Initial state of the tidal evolution of planetary systems 

The observed distribution of the synchronization parameter 
n/Q = Piot/PoTb provides us with information on the initial state 
of star-planet systems. We show n/Q vs. the eff'ective tempera- 
ture Teff of the host star for stars with Teff > 6000 K in Fig. [1] 
where the most impressive feature is the lack of systems in the 
range 1 < n/Q < 2, wi th the e xception of HAT-P-09 and XO-3. 

As conjectured bv iLanzal ([2010), this may be related to the 
processes occuring in the systems before their stars settle on the 
ZAMS. After the planet forms in the circumstellar disc, it mi- 
grates inwards by transferring its orbit al angular mom entum to 
the disc via resonant interactions (e.g., iLin et al.lll996l) . The mi- 
gration ends where the orbital period of the planet is half the 
orbital period of the particles at the inner edge of the disc where 
it is truncated by the magnetic field of the protostar. Since the 
protostellar magnetic field is also responsible for the locking of 
the stellar rotation to the inner edge of the disc, the value of 
the synchronization ratio at this stage is n/Q = PmtlPorh - 2. 
After a few million years, the disc disappears and the subsequent 
evolution depends on the interplay between tidal and magnetic 
interactions coupling the star and the planet. As the star con- 
tracts towards the ZAMS, it acc elerates its ro tation leading to 
a decrease of n/Q (cf. Fig. 8 of iLanza 1 l2010h . The loss of an- 
gular momentum through the stellar magnetized wind, however, 
counteracts such an acceleration. If the lifetime of the disc is 
i$ 10 - 15 Myr, the contraction prevails and the system enters 
into the domain n/Q < 2 characterized by a strong tidal interac- 
tion due to the dissipation of inertial waves in the extended con- 
vective envelope of the pre-main-sequence star (Ogilvie & Lin 
I2OO7I) . Angular momentum is transferred from the orbit to the 
spin of the star, leading to the infall of the planet towards the star 



(lHutlll980l) . Considering a star similar to the Sun, Q[ - 10 , a 
planet of Mo - 1 Mj, and an orbital period of 3 days, eq. (8) 
of Ogilvie & Lin (2007) gives an infall timescale of » 250 Myr, 
that may account for the lack of systems with 1 < n/Q < 2 given 
that typical stellar ages are between ~ 1 and ~ 8 Gyr. 

On the other hand, if the magnetic braking prevails, the sys- 
tem may be driven into the domain where n/Q > 2 and the tidal 
interaction becomes so weak {Q'^ ~ 10^-10'^) that the planet 
will spiral towards the star on a tim escale much loriger th an the 
main-sequence lifetime of the star ("Ogilvie. & Lin 1 120071) . thus 
populating the region of the diagram where n/Q > 2. Only 
when the planet is more massive than ~ 5 - 10 Mj and the star 
is late F or cooler, the tidal interaction may lead to the infall of 
the planet during the main-sequence lifetime of the star (cf., e.g., 
iBouchv et aOioilh . 

Another possible scenario occurs when the magnetic field of 
the pre-main-sequence star is so strong 3 - 5 kG) as to ef- 
fectively couple the orbital motion of the planet to the rotation 
of the star after the disc has disappeared dLovelace et al. 1120081 : 
IVidotto et al. Il201ft) . In this case, a synchronous state is attained 
on a timescale of ^ 2-3 Myr and maintained until the field 
strength decreases rapidly as the star gets close to the ZAMS. If 
the magnetic braking prevails over the final stellar contraction, 
the system would reach the ZAMS with 1 < n/Q < 2 and the 
planet will eventually fall into the star owing to the strong tidal 
interaction. On the contrary, if the stellar contraction prevails, 
the star will spin up, reach n/Q < 1 on the ZAMS, and the tides 
will push the planet away from the star. This may explain why 
several systems populate the domain n/Q < 1 in the diagram. 
The present value of n/Q is smaller than the initial value on the 
ZAMS because the tidal interaction transfers angular momentum 
from the stellar spin to the planetary orbit leading to an increase 
of the semimajor axis. Thus we shall base our calculations on 
the hypothesis that the systems observed today with n/Q < 1, as 
CoRoT-11, have reached the ZAMS in a state close to synchro- 
nization. 

Finally, remarks may be done concerning the two systems 
that have 1 < n/Q < 2. The peculiar state of the HAT-P-9 
system having n/Q ^ 1.4 may be explained by the weakness 
of the tidal interaction owing to the relatively small planetary 
mass (~ 0.7 Mj), and/or the young age of the system. For XO-3, 
the system has significant measured eccentricity and projected 
obliquity (e = 0.288 and A = 37.3 ± 3.7° respectively), so it is 
possible that the rotation period of the star is ill estimated or that 
the system has undergone a different evolution. 



3. Application to CoRoT-11 

3.1. System parameters 



i Gandolfi et al. I (1201 oh present the observations leading to the 
discovery of CoRoT-11 and derive the parameters of the system 
that are reported in their Table 4 together with their uncertain- 
ties. The stellar rotation period can be inferred only from the 
observed v sin / and the estimated radius of the star because no 
evidence of rotational modulation of the optical flux has been 
found in CoRoT photometry. As a matter of fact, a modulation 
with a period of » 8 days appears in the photometric time se- 
ries but it is incompatible with the maximum rotation period of 
1.73 ± 0.25 days as derived from the vsin/ = 40.0 + 5 km s"' 
and the star radius 7? = 1 .37 + 0.03 R©. It is likely to arise from a 
nearby contaminant star about 2 arcsec from CoRoT-1 1 that falls 
inside the CoRoT photometric mask. 
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The eccentricity of the planetary orbit is ill-constrained by 
the radial velocity curve due to the uncertainty of the radial ve- 
locity measurements for such a rapidly rotating F-type star. An 
upper limit of 0.2 comes from the modelling of the transit light 
curve that would lead to a mean stellar density incompatible with 
a F6V star for greater values. In any case, the present radial ve- 
locity and transit data are fully compatible with a circular orbit. 

The angle A between the projections of the stellar spin 
and the orbital angular momentum on the plane of the sky 
can be measured through the Rossiter-McLaughlin effect (e.g., 
lOhta et aTni2005h . Owing to the hmited p recision of the avail- 
able measurements, Gandolfi et al. '(2010) could not measure A, 
but they established that the planet orbits in a prograde direc- 
tion, i.e., in the same direction as the stellar rotation, and that 
A $ 50° - 80°. An aligned system, i.e., having A - 0°, is compat- 
ible with their observations. The obliquity e of the stellar equa- 
tor with respect to the orbital plane of the planet is given by: 
cos e = cos /p cos / + sin /p sin / cos A, where /p and / are the an- 
gles of inclination of the orbital plane and the stellar equator to 
the plane of the sky. Since /p = 83°3, cos e ^ sin / cos A, thus the 
measurement of A provides us with a lower limit for the obliquity 
of the stellar rotation axis. 



3.2. Stellar magnetic braking 



Late-type stars lose angular momentum thr ough their magne- 
tized stellar winds. In the case of F-type stars. lBarker & Qgilvie I 
( 1200 9) assumed a braking law of the form: ((iQ/(if)mb = -Tw^i"', 
where - 1-5 x 10"'^ yr is a factor that accounts for the effi- 
ciency of magnetic braking. It can be integrated to give: Q.{t) - 
fio[l + Twi^()(f - fo)] where Qo is the angular velocity at ini- 
tial time fo. For Q <K Qo this yields the usual Skumanich brak- 
ing law of stellar rotation. Applying that equation to CoRoT-lla 
leads to an unrealistically fast initial rotation for an age between 
I and 3 Gyr As a matter of fact, the above law of angular mo- 
mentum loss does not apply to fast rotators such as CoRoT-lla 
for which a saturation of the angular momentum loss rate is re- 
quired to expl ain the observ ed distribution of stellar rotation pe- 
riods (cf., e.g. iBouvier et al. il997 ). Therefore, we shall assume: 
(dQ./dt)m\, - -yw^^sat^' where Qsat is the angular velocity cor- 
responding to the transition to the unsaturated regime. The value 
of Qsat for solar-mass stars is » SQp, where Qp is the angular ve- 
locity of the present Sun (llrwin & Bouvier ' 2009). However, its 
value is highly uncertain in the case of mid-F stars because they 
do not show a rotational modulation measurable from the ground 
and experience a remarkably lower magnetic braking than G- 
type s tars. Considering the data in Table 4 of Wolff & Simon 
(Il997h . we estimate a rotation period corresponding to satura- 
tion of Si 2.8 days for mid-F stars with vsin/ ~ 20 - 25 km s"' 
at the age of the Hyades and the adopted value of yw However, 
if CoRoT-lla is older than ~ 2 Gyr, the value of the satura- 
tion period must be longer by a factor of at least » 2 otherwise 
the star would have been rotating faster than the brake-up veloc- 
ity on the ZAMS. This reduction of the angular momentum loss 
rate in the saturated regime can be a consequence of the initially 
fast angular velocity that led to a supersaturated dynamo regime 
for CoRoT-lla, or an effect of the close-in massive planet that 
reduced the efficiency of t he stellar w ind to extract angular mo- 
mentum from the star (cf. iLanza 1 12010) . 
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Fig. 2. The total angular momentum of the CoRoT-11 system in 
un its o f the cr itical angular momentum as defined in Eq. (4.1) 
oflHui (fT980l) versus the orbital mean motion in units of the 
present mean motion «obs- The stellar magnetic braking law has 
f sat = 6.0 days. The black solid line indicates the total angu- 
lar momentum values coiTesponding to tidal equilibrium. The 
black dotted line represe nt the stabil ity limit for the equilibrium 
as given by Eq. (4.2) of iHull ( Il98 0l). The solid blue lines show 
the backwards evolution of the system for different values of the 
tidal quality factor Q'^. From left to right, the values of Q'^ are: 
5 X 10\ 3 X 10\ 2 xlO\ 10\ 3.8 x 10^ and 10^ respectively. 
The initial points of the evolution are marked by the filled blue 
dots. 
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Fig. 3. The backward evolution of the orbital period f orb (solid 
line) and the rotation period (dotted line) vs. time since the 
present epoch for different values of the modified tidal quality 
factor Q'^ and a stellar braking law with f sat - 6.0 days. From 
left to right, the values of Q'^ are: 5 x 1Q\ 2 x IQ\ IQ\ 3.8 x 10^ 
and 10^, respectively. 



3.3. Backward integration from the present state to the initial 
condition 

We integrate the tidal evolution equations in Sect. 2.2 of 
iLeconte et al. backwards in time starting from the 

present condition, i.e., with the following parameters: stellar ro- 
tation period fiot - 1.73 days, planet orbital period Porb = 
2.99433 days, star mass M = 1.27 Mq, star radius R= 1.37 Rq, 
planet mass Mp = 2.33 Mj, planet radius Rp = 1.43 Rj (where 
Mj and Rj are the mass and radius of Jupiter), orbital eccentricity 
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Fig. 4. The backward evolution of the orbital period Poib (solid 
line) and the rotation period Prot (dotted line) vs. time since the 
present epoch for different values of the modified tidal quality 
factor 2s and a stellar braking law with P^at = 8.0 days. From 
left to right, the values of Q'^ ai-e: 2 x 10^ 1.4 x 10\ 8 x 10^ 
3 X 10^, and 10^, respectively. 



e = 0, and stellar obliquity e = 0°. The stellar quality factor Q'^ 
is regarded as a free parameter (see also Sect. 12.21 ). We account 
for the angular momentum loss of the star due to its magnetized 
wind by adopting a saturated loss rate with two possible values 
for f sat - 27r/Qsat, i.e., 6.0 or 8.0 days. Note that the results do 
not depend on the planet quality factor Q'^, provided that the or- 
bital eccentricity is zero and the planet is assumed to be in the 
synchronous rotation state. 

The present state of the system and its backwards evol ution 
are rep resented in a tidal stability diagram similar to that of lHutI 
(1198 0') in Fig.|2] The black solid line represents the value of the 
total angular momentum of the system when tidal equilibrium 
is established. For each value of the tot al angular momentum 
exceeding a minimum critical value (see lHutlll98d) . two equi- 
librium states are possible. The one on the left of the dotted line 
is stable, the other on the right is unstable. The filled square indi- 
cates the present state of the CoRoT-11 system with the assump- 
tion that ; = 90°, i.e., the stellar rotation period is 1.73 days. 

The results of the backwards integrations of the tidal equa- 
tions are represented by the solid blue lines for different values 
of Q'^ and a braking law with Psat - 6.0 days. All the backwards 
evolutions asymptotically approach an initial equilibrium state 
where all the time derivatives of the relevant variables vanish. 
Such an equilibrium is unstable, i.e., it leads to the transfer of an- 
gular momentum from the stellar spin to the planet that is pushed 
outwards. Increasing the value of makes the tidal interaction 
weaker so that the evolution of the total angular momentum is 
dominated by the magnetic braking of the star. This leads to an 
almost vertical evolution track in the diagram since the star was 
initially rotating faster There is no possible evolution track that 
crosses the dashed line, indicating that the past system evolution 
has occured only in the unstable region of the diagram, i.e., the 
planet has been pushed outwards starting from its initial orbit 
with «/Q < 1 during all the system lifetime. 

Given that the initial state is unstable, integrating the tidal 
evolution equations backwards in time is recommended to en- 
sure the stability of the obtained solutions. The instability of the 
present tidal state is not a specific property of CoR oT- 1 1 , but it is 
gene rally found in the case of hot Jupiter systems dLevrard et al. I 
l2009l) most of which evolve towards the infall of the planet into 



the star Therefore, what makes CoRoT-11 peculiar is the small 
value of the n/Q ratio pushing the planet outwards. 

The orbital period of the planet and the rotation period of the 
star are plotted vs. time in Fig.[3]for f sat - 6.0 days. If the initial 
state of the system was characterized by a synchronization of 
the stellar rotation with the planet orbital period, it is possible to 
constrain the tidal quality factor, i.e., 3.8 x 10^ < < 2 x 10^. 
Although we speculated that the system may have reached the 
ZAMS close to synchronization in Sect. 12.41 other evolutionary 
scenaria are possible. In some of these scenaria, the planet may 
have started with an orbital period longer than the stellar rotation 
period 3 Gyr ago, leading to the present state if > 2 x 10^. 
As an example of such an evolution, we plot in Fig.|3]the case of 
Q's - 5x 10^. Note that in this case the initial state of the system 
has n/Q. ^ 0.25 that is significanly smaller than the observed 
values of n/Q. plotted in Fig. [1] Although the present sample 
of systems with n/Q < 1 is too small to draw conclusions on 
the minimum value of n/Q during their evolution, we hope that 
such a statistical information can be used to constrain Q'^ when 
a larger sample is available. 

On the other hand, Q'^ < 3.8 x 10'' implies a system younger 
than ~ 1 Gyr. To illustrate this point, we plot in Fig.[3]the evolu- 
tion for 2s - 10^ for which the system reached synchronization 
~ 0.4 Gyr ago. Although the backwards evolution can be con- 
tinued to earlier times, this initial state is unstable as well as all 
the allowed initial states with Porb > ^'rot- For Ql < 3.8 x 10^, 
the tidal interaction would be so strong and the evolution so fast 
that the system would reach a present value of Poib longer than 
observed. Therefore, we can set a lower limit on Ql on the ba- 
sis of the estimated age of the system. Note that if the system 
were younger than ~ 0.7 - 1 Gyr, the Q'^ value could be so small 
that we may directly measure it after a few decades of timing 
observations (see Sect.|4]i. 

The evolution tracks obtained for f sat = 8.0 days and differ- 
ent values of the average Q'^ are plotted in Fig.|4] They show that 
a decrease of the braking rate of CoRoT-lla by ~ 30 percent 
leads to a comparable decrease in the values of the average Ql 
because a stronger tidal interaction is needed to transfer the re- 
quired amount of angular momentum from the stellar spin to the 
planet when the stellar rotation was initially slower. On the other 
hand, if we adopt Psat - 2.8 days, as estimated from the data of 
Wolff & Simon (1997), the stellar braking was so rapid that only 
in the case of a system younger than 1 Gyr we may find an ac- 
ceptable evolution with Q'^ ~ 9 x lO*". For higher values of Q'^, 
the evolution takes longer and this leads to a planet approaching 
the Roche lobe radius and/or an initial stellar rotation exceed- 
ing the break-up velocity. Therefore, the assumption that Psat is 
at least ~ 5 - 6 days seems to be more likely for CoRoT-lla 
than Psat ~ 3 - 4 days. In the latter case, the system would be 
so young and the average Q'^ so small that we expect to detect 
the tidal increase of the orbital period within a few decades (see 
Sect. nil. 

3.3.1 . Initial eccentricity of the orbit 

In the scenario introduced in Sect. 12.41 the initial eccentricity 
of the planetary orbit would be close to zero ( Sari & GoldreichJ 
2004), unless an encounter with some planetary or stellar body 
would have excited it. Assuming that no third body has been 
perturbing the system after that event, the evolution of any 
primordial eccentricity is ruled by the dissipation inside both 
the star and the planet, the latter being the most important for 
Q'j, ^ 10**. The timescale for the damping of the eccentricity 
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Fig. 5. The forward evolution of the initial obliquity from the 
ZAMS for different initial value of the angle itself. The planetary 
orbit is assumed circular and the initial orbital period is 0.7 days, 
while the initial rotation period of the star is 0.48 days for e(f = 
0) = 45°, 0.60 days for e{t = 0) = 30°, and 0.67 days for e(f = 
0) = 15°, to prevent the infall of the planet. 



Tg = { \ I e){de I dt) adopting the present parameters of the sys- 
tem is ~ 120 Myr for e = 0.2, = 5 x 10^ and Q'^ = 10^ 

and increases to ~ 1.1 Gyr for Q'^ - 10^. Since the tidal interac- 
tion was initially stronger - it scales as a^^ - the initial Tg was 
one order of magnitude shorter than the above values, leading 
to a rapid damping of any initial eccentricity, unless Q'^ ^ 10^. 

Moreover, for initially large eccentricities (i.e., e ^ 0.3) and an 
initial n/Q. x 1, the evolution of the system may be significantly 
affected and the infall of the planet occurs as the orbital angular 
velocity at the periastron exceeds the stellar angular velocity. To 
prevent such an infall, we must assume that the star was initially 
rotating remarkably faster than considered above. Since the ro- 
tation braking on the main sequence is expected to be small, this 
is hardly compatible with the present estimated rotation period. 
We conclude, that the present eccentricity of CoRoT-1 1 is likely 
to be close to zero. In any case, in Sect. |4] we shall describe ob- 
servations that can constrain the present value of e to allow us to 
discriminate among different scenarios for the system evolution. 

3.3.2. Initial obliquity of the orbit 

The obliquity of the CoRoT-1 1 system may be significantly dif- 
ferent from zero, as indeed found for seve ral stars accompanied 
by hot Jupiters with T^ff > 6250 K ( Winn et al. II2OIO) . In the 
case of a system born close to n/Q = 1 , we need Q cos e > n to 
prevent the tidal infall of the planet and the observational result 
that CoRoT-1 1 is compatible with an aligned system is in agree- 
ment with this constraint. Initial values of the obliquity e > 30° 
require a fast-rotating star which is difficult to reconcile with the 
presently observed vsin/ and stellar radius. However, smaller 
values of the obliquity cannot be excluded. They are damped 
on a timescale that depends on the value of Q'^, as illustrated 
in Fig. |5] where the time has been scaled inversely to the stel- 
lar tidal quality factor. These plots show that, after an initial re- 
markable damping of the obliquity, it asymptotically approaches 
a value about half its initial value that pratically corresponds to 
the present value. In Sect. |4]we shall discuss how such a value 
can be constrained with accurate timing observations of the sys- 
tem. 



4. Measuring the present tidal dissipation 

The main conclusion of the previous section is that the mean 
modified tidal quahty factor of CoRoT-lla < 2 x 10^, pro- 
vided that the system started its evolution on the ZAMS from 
a nearly synchronous state. Since Q'^ is subject to remarkable 
variations during the evolution of the system, this result is of 
limited value to constrain dynamical tide theory, as in the case 
of the mean values of Ql derived from close binaries in clus- 
ters. However, a much more interesting opportunity offered by 
CoRoT-11 is the possibility to measure the present value of Q'^ 
from the rate of orbital period variation. Moreover, as we shall 
see, the rate of orbital period variation depends also on the or- 
bit parameters, and we can extract a wealth of information from 
CoRoT-1 1 using the timing method. 

4.1. The circularized and aligned system hypothesis 

If the present orbit is circular, the semimajor axis is subject to a 
secular increase because the star is rotating faster than the orbital 
motion of the planet, so angular momentum is transferred from 
the stellar spin to the orbit via the tidal interaction. When e = 
and the stellar obliquity is negligible (e ~ 0°), the rate of such 
a transfer depends on the system parameters, that are known to 
a good level of precision, and on Q'^ only. Therefore, we can 
compute the rate of orbital period increase dPoih/dt and use it 
to derive the difference {O - C) between the observed time of 
mid-transit and that computed with a constant period ephemeris 
after orbital periods, viz.: 

(O - C)(N) = i(^)P,^N(N + 1), (3) 
2 dt 

where we have assumed that the constant period of the 
ephemeris is the orbital period at A^ = 0, i.e., at t he epoch of the 
first tr ansit. Using the tidal evolution model of iLeconte et al. I 
and assuming e = 0, e = 0°, and the maximum stellar 
rotation period Piot = 1 .73 days, we find: 

(0-0 = 8.34(^1' s (4) 

after N - 3000 orbital periods, i.e., 24.6 years. 

The epoch of mid transit can be determined with an accuracy 
that depends on the depth and duration of the transit, the accu- 
racy of the photometry, and the presence of noise sources, e.g., 
related to stellar activity. In the case of CoRoT-1 1, since the star 
shows no signs of magnetic activity, the precisi on of the transit 
timing is dominated by the photometric accuracv. lGandolfi et al. I 
(2010) report in their Table 4 a precision of ~ 25 s for the initial 
epoch of their ephemeris as obtained from a sequence of ~ 50 
transits observed with a photometric accuracy of ~ 250 p arts per 
million (ppm) and a cadence of ~ 130 s. Recently, Win n et aTl 
(2009) using the Magellan/Baade 6.5-m telescope, observed two 
transits of the hot Jupiter WASP-4b, reaching a photometric ac- 
curacy of ~ 500 ppm (standard deviation) in their individual 
measurements obtained with a cadence of 30 s. The average stan- 
dard deviation of the epoch of mid-transit is ~ 5 s which can 
be assumed as the limit presently attainable with ground-based 
photometry. The transit of WASP-4b has a depth of ^ 0.028 mag 
and the star has a magnitude of V ^ 12.5, while that of CoRoT- 
11 has a depth of x 0.011 mag and the star has V = 12.94. 
Therefore, it is reasonable to assume that a timing accuracy of 
0.5 -Is can be reached with a space-borne dedicated telescope 
of aperture 2 - 3 m, since the precision is expected to scale with 
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the accuracy of the photometry that should reach 200 - 300 ppm 
in ~ 15 s integration time. Such an instrument would allow us to 
detect the expected transit time variation (hereafter TTV) with 
3cr significance after an interval of ~ 25 years for Ql ~ 5 x 10^. 
Even in the case that no significant TTV is detected, a stringent 
lower limit can be placed on Q'^, with the sensitivity increasing 
as the square of the considered time baseline (cf. Eq.|3]l. 

4.2. The eccentric system hypothesis 

A significant eccentricity of the planetary orbit can affect the 
O - C of the TTV, particularly when tidal dissipation in the 
planet is not negligible. This happens because the angular mo- 
mentum of the planetary orbit is given by Lp - [(MsMp/(Ms + 

so an increase of Lp does not necessarily lead 
to an increase of a when the eccentricity is damped. Assuming 
e = 0.2, = 5 X 10^ and = 10^ we find O - C = -8.33 s 
after 24.6 yr, owing to the strong dissipation inside the planet 
that dumps the eccentricity with a timescale of ~ 120 Myr. If 
the dissipation inside the planet is smaller, the absolute value 
of the O - C decreases. For instance, if Q'^ = 10^, we obtain 
O - C = 0.49 s after 24.6 yr 

Furthermore, if the orbit of CoRoT-11 is indeed eccentric, 
the measurement of the O - C of the primary transit are not 
sufficient to determine Q'^. This happens because in addition to 
the tidal evolution, also the precession of the line of the apsides 
(due to the planetary and stellar tidal de formations and general 
relativistic effects) pro duces a TTV (see lMiralda-Escudell2002t 
iRagozzine & Wolf 2009). For e - 0.2, and a planet Love number 
k2 - 0.3, similar to that of Saturn, we have a precession period 
of only ~ 3300 yr, that leads to a maximum (9 - C = 18.1 s after 
A' = 3000 orbital periods, i.e., 24.6 yr. To separate the effects 
of the apsidal precession from those of tidal origin, we note that 
the apsidal precession produces an oscillation of the 0-C of the 
mid of the secondary eclipse, i.e., when the planet is occulted by 
the star, that is always shifted in phase by 180° with respect to 
that of the primary transit. In other words, the two oscillations 
have always opposite signs as it is well known i n the case of the 
apsid al precession of close binaries (see, e.g., iGimenez et al. I 
Il987h . On the other hand, the O - C of tidal origin has the same 
sign for both the primary transit and the secondary eclipse, so it 
can be separated from that arising from apsidal precession. 

We conclude that a measurement of the eccentricity (or a de- 
termination of a more stringent upper limit) is critical to derive 

from the 0-C observations. A direct constraint on the ec- 
centricity can be obtai ned from the o bservation of the time of 
the secondary eclipse. Ga ndolfi et al. I (12010) did not find evi- 
dence of a secondary eclipse in the CoRoT light curve which 
sets an upper limit of ~ 100 ppm to its depth in the spectral range 
300 - 1 100 nm. However, it should be possible to detect the sec- 
ondary eclipse in the infrared, thus deriving the value of e cos w, 
where o) is the longitude of the periastron, from the deviation of 
the time of the secondary eclipse from half the orbital period. In 
the case of HD 189733, ob serving in a band centr ed at 8 yum with 
IRAC onboard of Spitzer, Knutson e t al. I (12007 ) found that the 
secondary eclipse occurred 120 ± 24 s after the time expected for 
a circular orbit, leading to e cos co - 0.001 ± 0.0002. The same 
method has been applied to constrain the eccentricit y of the orbit 
of CoRoT-2 finding e cos w = -0.0029 1 + 0.00063 ( Gillon et al. I 
12010 ). Moreover, the duration of the primary transit and the sec- 
ondary eclipse can provide a measurement of e sin o) although 
with a precision significantly smaller than in the case of the sec - 
ondary eclipse timing (see Sect. 2.5 of lRagozzine & Woill2009h . 



In conclusion, we recommend to perform infrared photometry of 
CoRoT-1 1 to constrain its orbital eccentricity. 

We note that the present fast rotation of CoRoT-lla can in 
principle excite the eccentricity of the orbit during the phase 
of the evolution w hen Prot/Porh < 11/18 (cf. Eq. (6) of 
iLeconte et al. ll20Toh . However, the tidal dissipation inside the 
planet is effective at damping the eccentricity and only for Q'^ > 
2.5 X 10^ we find an increasing eccentricity for the present state 
of the system. Imposing that the eccentricity is always damped 
along the past evolution of the system, we find Qp < 1.6 x 10^ 

for > 4 X lO''. 

4.3. Non-negligible obliquity hypothesis 

If the present obliquity of the system is not negligible, the rate 
of orbital period increase and the corresponding 0-C in Eq. (|4|i 
must be multiplied by: 



because the component of the stellar angular momentum inter- 
acting with the orbit is cos e, where Lj is the modulus of the 
angular momentum of the star (cf. Eq. (2) of iLeconte et al. I 
2010). 

It is interesting to note that the obliquity can be constrained 
not only through the observation of the Rossiter-McLaughlin ef- 
fect, but also by measuring the precession of the orbital plane of 
the planet around the total angular momentum of the system that 
is expected when e 9^ 0°. Such a precession arises because the 
stellar gravitational quadrupole moment is non-negligible due to 
its fast rotation. M ore precisely, adopting a stellar Love number 
k2s = 0.0138 after lClaretl (Il995h and P^ot = 1-73 days, we find 
J2 = 4.19 X 10"^. The contribution to J2 arising from the tidal 
deformation of the star induced by the planet is about 1 100 times 
smaller and can be safely neglected. 

Even assuming an obliquity of a few degrees, the precession 
period is so short as to give measurable effects. For instance, 
for e = 5°, the precession of the orbital plane has a period of 
~ 7.0x10"^ years, according to Eq. (4) of Miralda-Escude (2002?). 
This produces a variation of the length of the transit chord across 
the disc of the star leading to a variation of the transit duration 
(transit duration variati on, hereafter TDV). Applying Eq. (12) of 
Miralda- Escudi(l2002h and assuming that the total angular mo- 
mentum of the system is approximately perpendicular to the line 
of sight, we find a maximum TDV of ~ 8 s in one year. Assuming 
a precision of ~ 40 s in the measurement of the duration of the 
transit as obtained with CoRoT observations, such a variation 
should be easily detected in ~ 5 yr and can provide a lower limit 
for the stellar obliquity, remarkably better than achievable with 
the Rossiter-McLaughUn effect if e < 15° - 20°. 

5. Role of a possible third body 

The use of TTV to measure the tidal dissipation inside the star 
(and the planet, in the case of an eccentric orbit), can be ham- 
pered by the presence of a distant third body that induces appar- 
ent 0-C variations due to the light-time effect related to the 
revolution of the star-planet system around the common centre 
of mass. In principle, any light-time effect induced by a third 
body is periodic while the effects of tidal dissipation are secu- 
lar, i.e., they do not change sign. However, if the orbital period 
of the third body around the barycentre of the system is longer 
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than, say, three times the time basehne of our observations, it 
may become difiicuh or impossible to separate the two effects. 

To be specific, let us consider an (9 - C of 3 s observed 
with a time baseline of AT ~ 25 yr. If it is due to a light- 
time effect, the corresponding radial motion of the star-planet 
system is Az = c{0 - C) ~ 10^ m, where c is the speed of 
light, and the radial velocity of the barycentre of the system 
is Az/AT ^ 1.26 m s"', too small t o be detectable with the 
present radial velocity accuracy (cf. iGandolfi et al. ll20Toh . A 
third body moving on a coplanar circular orbit with a period 
of, say, Ptb = 75 yr, would have an orbital semimajor axis of 
~ 2.88 X 10'^ m and a minimum mass of 0.2 Mj to induce a 
reflex motion of the above amplitude. Its direct detection is very 
difficult since its apparent separation from the central star would 
be only 0.034 arcsec, given an estimated distance to CoRoT-11 
of 560±30 pc. Astrometric measurements with the GAIA or SIM 
satellites cannot be used to detect the reflex motion of CoRoT-1 1 
since it would correspond to a maximum displacement of only 
~ 2.4 yuarcsec during an interval of 5 yr, that is impossible to 
separate from the much larger proper motion of the system. 

If the orbit of the third body is inclined to the orbital plane of 
the star-planet system, the orbital plane of the planet will precede 
around the total angular momentum of the system. Assuming a 
third body mass of 0.2 Mj, a relative inclination of the orbits of 
45°, and a semimajor a xis of ~ 2. 88 x 10'^ m, equal to ~ 20 AU, 
Eq. (8) of Miralda-EsciH (12002 ) gives a precession period of 
~ 14 Gyr, i.e., too long to be observable. 

Secular perturbations of the orbit of the planet by the distant 
third body can be neglected because they induce a variation of 
the O - C of the order of Pprb/f tb ~ 10 of the light-time eff'ect 
(cf., e.g., iBorkovits et al. Il2003h . On the other hand, a perturber 
on a close-in orbit, possibly in resonance with the planet, would 
induce O - C variations with a short periodicity and an ampli- 
tude up to a few minutes ( Agol et al. 2005), that would be easily 
detectable. 

It is interesting to note that CoRoT- 11a has a low level of 
magnetic activity in contrast to the case of stars with lower effec- 
tive temperature. This implies that the orbital period modulation 
induced by magnetic activity (Lanza et al, 1998; Lanza 2006) 
can be excluded as a source of 0-C variations in this c ase. This 
is not possible for cooler stars (* Watson & Marsh Il20 10l) that are 
therefore less suitable to measure the TTV induced by tidal dis- 
sipation. 

6. Conclusions 

Assuming an initial state of the CoRoT-11 system close to syn- 
chronization between the stellar spin and the orbital period of 
the planet, we can put constraints on the average modified tidal 
quality factor of its F-type star, finding 4xl0^$2s~2x 
10^. Rigorously speaking, we should have expressed these con- 
straints in terms of the average tidal time lag At because our 
tidal equations are valid for a constant At (cf. Sect. |2] and 
iLeconte et al. Il20l"ol) . However, in view of the uncertainty on the 
above constraints, we prefer to give them in terms of that 
varies by a factor up to a; 4-5 during the evolution, if Af is 
assumed to be constant (cf. Sects. I2land l3. 3b . 

An initial state close to synchronization is not the only pos- 
sible one. However, the minimum estimated age of the system 
of ~ 1 Gyr implies Ql^4 x 10^, otherwise the evolution from 
any initial state with Porb > Aot would be too fast, leading to a 
present Poih longer than observed. An upper limit on Q'^ can be 
set when a larger statistics of the values of the n/Q. - PmtlPorh 



ratio will be available because for 2^ ^ 5 x 10^ the initial value 
of n/Q is predicted to be ~ 0.25 that is smaller than the presently 
observed minimum n/Q. ~ 0.5 (see Sect. l3.3l l. 

The above limits on Q'^ are only a factor of 2 - 15 s maller 
than the average value predicted by iBarker & Oplvie I (l2009h 
which can be consider ed a success of the dynamic tide theory of 
Qgilvie & Linl (|200"A given the current uncertainty on the tidal 
dissipation processes occurring in stellar convection zones. 

We find that CoRoT-1 1 is one of the best candidates to look 
for orbital period variations related to tidal evolution by moni- 
toring its transits and secondary eclipses in the optical and in the 
infrared passbands. In contrast to several systems hosting very 
hot Jupiters, e.g., WASP-12 or W ASP- 18, whose orbital decay 
may be too slow to be measurable vie &Linl2007tlLietal.l 
2010), the parameters of CoRoT-1 1 appear to be ideal for a de- 
tection by measuring the times of mid-transit and mid-eclipse 
with an accuracy of at least 1 - 5 s along a time baseline of 
~ 25 yr This is a consequence of the peculiar ratio of the ro- 
tation period of the star to the orbital period of the planet that 
is presently the smallest among systems having a radial velocity 
curve compatible with zero eccentricity. The test is particularly 
sensitive to values of Q'^ between 10^ and lO*" leading to a fast 
tidal evolution of the system, as discussed in Sect. 13.31 

The CoRoT-1 1 system is also particular because the preces- 
sion of the orbital plane due to a non-zero obliquity of the stel- 
lar spin can be measured through the variation of the duration 
of the transit on a timescale as short as 5 - 10 yr by means of 
ground-based observations with telescopes of the 6 - 8 m class. 
Note that even an obliquity as small as 5° can be detected with 
this method. On the other hand, the eccentricity of the planetary 
orbit can be well constrained by measuring the times of the sec- 
ondary eclipses in the infrared and their duration. If the system 
turns out to have an eccentric orbit, this puts a constraint also 
on the minimum value of the planet quality factor Q'^ if we as- 
sume that the eccentricity is of primordial origin, or is a clear 
indic ation of the presence o f a perturbing body that excites it 
(e.g., iTakeda & Rasio 1 l2005h . The light-time effect due to a dis- 
tant third body is indeed the only phenomenon that can seriously 
hamper the detection of the transit time variations expected from 
the tidal orbital evolution. The detection of a distant companion 
may be particularly difficult, but this limitation is present also in 
any other system candidate for a direct measurement of the tidal 
dissipation. 
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